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Quantum critical systems out of equilibrium are of extensive interest, but are difficult to study the- 
oretically. We consider here the steady state limit of a single-electron transistor with ferromagnetic 
leads. In equilibrium (i.e., bias voltage V = 0), this system features a continuous quantum phase 
transition with a critical destruction of the Kondo effect. We construct an exact quantum Boltz- 
mann treatment in a dynamical large-N limit, and determine the universal scaling functions of both 
the non-linear conductance and fluctuation-dissipation ratios. We also elucidate the decoherence 
properties as encoded in the local spin response. 
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Quantum phase transitions are of extensive current in- 
terest in a variety of strongly correlated electronic and 
atomic systems A quantum critical point occurs when 
such transitions are second order Q. In a quantum crit- 
ical state, there is no intrinsic energy scale in the exci- 
tation spectrum. External probes will readily drive the 
system out of equilibrium, such that the standard notion 
of linear response breaks down. The response of such a 
system to an external drive cannot be directly related to 
its intrinsic fluctuations in equilibrium. 

Compared to their classical counterparts, quantum 
critical systems are more difficult to study already at the 
equilibrium level. A quantum critical point involves the 
mixing of statics and dynamics. This complicates the de- 
termination of the equilibrium fluctuation spectrum, es- 
pecially in the long-time "quantum relaxational" regime 
(thus <C kgT) at finite temperatures 0. In addition, the 
classification of the universality classes is delicate, since 
the critical modes may incorporate inherently quantum 
ones that go beyond the order-parameter fluctuations. A 
typical example is for Kondo systems, in which quantum 
criticality can involve a critical destruction of the Kondo 
entanglement 0. 

To make progress theoretically, we need approaches 
that can, on the one hand, capture the scaling prop- 
erties of a quantum critical point in equilibrium, and, 
on the other hand, be extended to settings out of equi- 
librium. This is impeded by the lack of a general un- 
derstanding of how to generalize the free energy func- 
tional to stationary nonequilibrium states from which 
scaling relations could be obtained. Some progress 
has been made in recent years. In the case of the 
two-dimensional superconductor-insulator transition, the 
steady-state conductance in the non-linear regime [V^(s> 
ksT / e) — > 0] is found to differ from its counterpart in 
the linear-response regime \T{~^> eV/kB) — ► 0] [I, HI, in 
analogy to the non-commutativity of the us — > and 
T — > limits of the equilibrium current-current corre- 
lation function^. A complementary issue, viz. quan- 
tum criticality induced by a steady-state voltage drive, 



has also been addressed in the contexts of itinerant 
magnetism 0, 0|. Still, there is a considerable need for 
insights into many open issues regarding quantum criti- 
cal systems out of equilibrium. How do we connect the 
scaling behavior of the equilibrium fluctuation spectrum 
with that of the out-of-equilibrium properties? How to 
calculate the universal out-of-equilibrium scaling func- 
tions? Is the fluctuation-dissipation theorem violated in 
a universal way in the quantum critical regime? 

In this letter, we study these issues in the steady 
state of a quantum critical system out of equilibrium. 
We have chosen a ferromagnetic single-electron transis- 
tor (SET)0 as a case study. It was shown earlier that 
this system undergoes a gate-voltage-induced quantum 
phase transition, which is characterized by a continuous 
suppression of the Kondo screening [l2l [l3j . A dynamical 
large-N limit was found to capture the universal 
erties of the quantum criticality in cquilibriumf 
Here, we show that this approach can be reliably gen- 
eralized to the case of steady states under a finite bias, 
and use it to determine the full scaling functions of non- 
equilibrium properties. Our results establish the system 
as a prototype model setting in which the general issues 
raised earlier can be systematically studied. In addition 
to elucidating the nonequilibrium quantum criticality, 
our work also contributes to the general understanding of 
the Kondo effect out of equilibrium [II, [ll, [13, H, 0, lS]. 

The system, the model, and its large-N generaliza- 
tion: The system we consider is illustrated in Fig. la, 
showing a quantum dot coupled to two leads of fer- 
romagnetic metals whose magnetizations are antiparal- 
lel. It was demonstrated experimentally[9( and theo- 



retically [2l|, l22| that a Kondo effect can still occur 
in the Coulomb-blockade regime, in spite of a large 
Stoner splitting of the conduction electron bands. It was 
shown[12l,ll3| that the system should be modeled in terms 
of a Bose- Fermi Kondo Hamiltonian. A localized moment 
on the dot is coupled to a fermionic bath, corresponding 
to the Stoner continuum of the ferromagnetic leads, with 
a coupling constant J. Simultaneously, it is coupled to a 



bosonic bath, associated with the magnons of the leads, 
with a coupling constant g m J 2 /T (where T is the bare 
resonance width). Tuning the gate voltage varies the ra- 
tio (g/J) of the two couplings, giving rise to a continuous 
transition between the Kondo phase and a critical local 
moment (LM) phase. 

Of interest to us here is the universal behavior of the 
quantum critical point, referred to henceforth as C, and 
the quantum critical LM phase. The equilibrium behav- 
ior is illustrated in Fig. lb, for a fixed J (and, equiva- 
lently, a fixed bare Kondo scale Tj^); the quantum critical 
point C corresponds to g — g c , and the quantum critical 
LM phase occurs at g > g c . Universal scaling behavior in 
this model sets in whenever the probing energy and the 
temperature are smaller than min[T^-,A s ], where A s is 
the energy at which the magnon dispersion merges with 
the Stoner continuum. Physically, we have <c A s and 
we therefore restrict ourselves to T <C T^. 

Anticipating the large-N limit to be taken, we write 
the effective low-energy model in the following form: 

H = (J/N)S a ■ c\ a r a aa ,c la , + H Q (c) 



(1) 



Here, c[ a = (t L c\ aL + t R 4 aR )/ y/\t L \ 2 + \t R \ 2 creates 
an electron in a linear combination of the local states 
in the leads for spin a and channel index i ; a, a 1 = 
1, . . . , N, i = 1, . . . , M, and th and t R are the hybridiza- 
tion matrix elements with the left and right leads, re- 
spectively. Hq(c) = J2 a =L,R H 0a(Ca), with H 0a (c a ) = 

Efc«r( e fc,^ ~ /OcjL, a c fcio . )Q1 . The energy dispersion, 
£k,a,a, contains the appropriate Zeeman shifting, but is 
otherwise taken to be featureless and has in particular a 
finite density of states in between the two chemical po- 
tentials and \i R . 3? = Y^qi^l + 1S a ^ so a linear 
combination of the spin waves on the left and right leads; 
it has been written to contain N 2 — 1 components. The 
matrices r a ,a = 1, . . . , A^ 2 — 1 span the generators of 
SU(N). For simplicity, we consider t^ = t R (otherwise, 
a local magnetic field will arise [3, Hi]]). It follows from 
uj q cx g 2 that the spectrum of the bosonic modes is 

P(w - w g ) - S(u + w,)] ~ M^sgnM, (2) 



with e = 1/2. (Our results can be easily generalized to 
the generic sub-Ohmic case, < e < 1.) The power-law 
form given in Eq. @ is valid for lu below some cutoff fre- 
quency, A. For simplicity, we consider the particle-hole 
symmetric case without any potential-scattering term. 
The latter is exactly marg inalO EEI ~ 



Il2j . and would 
hence only add a constant offset without modifying the 
T- and ^-dependent scaling terms. The explicit Keldysh 
calculation in the presence of the potential scattering 
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FIG. 1: (a) A sketch of the magnetic SET. (b) The phase 
diagram of the magnetic transistor in equilibrium (fiR = /il). 
Tfc is the bare Kondo scale, appearing when g = 0, i.e., when 
there is only the coupling J. A critical point at coupling 
g c separates the Kondo phase from a critical local moment 
(LM) phase, (c) The Keldysh contour appropriate for the 
non-equilibrium cases. 



term is more difficult, because solving the pscudofermion 
constraint in the absence of particle-hole symmetry is 
more involved; we have nonetheless used such a calcula- 
tion to confirm the above[23j|. 

The dynamical large-N limit in equilibrium was consid- 
ered in [14j |. and this limit has been shown to correctly 
capture the dynamical scaling properties (including, in 
particular, the finite-temperature relaxational proper- 
ties) of the N = 2 case for both the quantum critical 
point (g = g c ) and quantum critical LM phase (g > 



5c) PJ, lij- In the large-N limit [2J, [25(, S is expressed 
in terms of pseudofermions, S aa i = f^f a i —S aa >q . where 
<7o = Q /N is related to the chosen irreducible representa- 
tion of SU(N) and it also appears in the pseudo-fermion 
occupancy constraint ~^2 a f^f a = Q. The quartic inter- 
action is decoupled via a Hubbard-Stratonovich field Bi 
(i = 1, . . . , M), which is conjugate to £ CT c Lfa/^- 

For a finite bias voltage, we work on the Keldysh con- 
tour, depicted in Fig. Since we are interested in the 
steady state limit, we specify the state of the system at 
the infinite past to = — oo, so that at finite time all initial 
correlations will have washed out. A finite bias voltage, 
eV = /il — Mflj is incorporated in the following lesser 
function for the fcrmionic bathlla. 12611 : 



Gf(fj) = m[f R (u) + f L {uj)]p(u), 



(3) 



with f L / R (x) = [exp((z - fi L/R )/T) + l] and ^ L = 
eV/2 = —p,R. Eq. ([3]) also defines the lead distribution 
function. Both leads are held at temperature T which 
completely characterizes the magnon distribution func- 
tion. The magnons, therefore, retain thermal equilib- 
rium, and the corresponding lesser function is 



G£ = -2irib(w)A&(u). 



(4) 



Here, = [l/r(l/2)]|w| 1 / 2 sgn(w), where T is the 

Euler Gamma function, and b(ui) = [exp(ui/T) — 1] 1 . 

Large-N limit on the Keldysh contour: In the equi- 
librium case, the large-N limit yields saddle-point equa- 
tions for the pseudoparticle self-energies 14{. On the 
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FIG. 2: Scaling of the current-voltage ('IV') characteristics 
near (a) the critical (C) fixed point (g = g c ) and (b) the 
critical LM phase (g = 5g c > g c ). In (a), we have shown the 
results over a reduced range of V/T in order to highlight the 
crossover in the scaling function; scaling is in fact seen over a 
similarly extended range (nearly 10 decades) of V/T as shown 
in (b). The results are shown with h,e,kB being set to 1. 



Keldysh contour, the saddle point equations become 
E / (t) = - K ig c (t)G B {t) - g 2 ig$(t)G f {t), and E B (i) = 
iQ c (t)Gf(—t), with k = M/N. Performing an analytical 
continuation onto the real time axisl27l 12811. we have 



E5(t) 



E>(i) 



-«»0<(*)G£(t) 
-K»$>(t)G£(t) 
ig?(t)G>(-t), 
iG>{t)Gf{-t). 



g 2 iQt{t)G<{t), 
9 2 iG£(t)G>(t), 



(5) 



Retarded (r) and advanced (a) functions are related to 
the lesser and greater functions via G a (t,t') = Q(t' — 
t)(G < (t'-t)-G > (t'-t)) and G r (t,t') = Q(t-t')(G> (t' - 
t) - G<(t' - t)), so that G r - G a = -2iIm(G) = G> - 
G<. G< = G r E<G a , where G r /G a are obtained from 
G-^lu) = G \uj) - E(w), with G-g = 1/J and G~j = 
w — A; A is a Lagrangian parameter enforcing the pscudo- 
fermion occupancy constraint. 

We solve the saddle point equations numerically for 
given temperature T and voltage V on real frequencies. 
We choose k = 1/2 and po(uj) = p cxp(— uj 2 /tt). The 
effective Kondo coupling is taken to be Jp = (Jll + 
Jrr)Po — 0.8/7T, resulting in a bare nominal Kondo scale 
Tftpo = 0.06/7T at g = 0. The bosonic bath is cut off 
smoothly at Apo = 0.05/7T. This yields a critical coupling 
g c po ph 0.79/7T. For the LM regime we will illustrate the 
result at g = 5<? c , which is sufficiently large to ensure 
an extensive scaling regime. The convergence criterion is 
that, for both G^ — Cj and G^ — G R , the sum over the 
entire frequency range of the relative difference between 
two consecutive iterations is less than 10 ~ 5 . 

Steady-state current: The current between the dot and 
the left lead for an arbitrary bias voltage V is [29[ 



h = - J dup(u,)[f L ( U )(7Z L ( U )-7Z L {v))+T< L (w)] 
In the steady state limit, II + Ir = 0, and the current 



can then be alternatively cast into: 

I(V,T) = ^ J <kjp(u)[f L (u)-f R (u)} 

x 1m[T£ L (u>,T,V) +T£ R (u>,T,V)]. (6) 

Here, T a ,p, («,/? = L/R), is the T-matrix of the Bosc- 
Fermi Kondo model. In the large-N limit we consider, 
T{t) = T LL + T RR = (i/N) G B (-t)G f (t). 

V/T scaling of the I-V characteristics: We are now 
in a position to discuss the quantum critical behavior 
out of equilibrium. Fig. [D^a) shows the current-voltage 
characteristics at the quantum critical point, g = g c . 
In the linear-response regime, for T — > after taking 
V — > 0, the conductance is a universal constant in the 
leading order[I3|. In the non-linear regime, for V — ► at 
T = 0, we find that the conductance is also a universal 
constant, but differs from the linear-response value. This 
is reminiscent of what happens Q at the two-dimensional 
superconductor-insulator transition. It is to be con- 
trasted with the result in the two-channel Kondo prob- 
lem, in which the two constants are identical; V/T scal- 
ing is then analyzed in an unorthodox fashion, in terms 
of AG = G(V,T) - G(V = 0,T = 0) 0, SH Hi . 



Beyond these limits, we are able to determine the con- 
ductance for arbitrary V, T. We find that I/T shows a 
scaling collapse of V/T, provided V,T T£. 

Consider next the quantum critical LM phase at g > 
g c . In the linear response regime, the conductance is 
proportional to T 1 / 2 12j. In the non-linear regime, V — > 
at T — 0, we find that the conductance goes as V 1 / 2 . 
Fig. now demonstrates a scaling collapse of I/T 3 / 2 
vs. V/T. For both cases, Fig. 2 shows the full scaling 
function for the non-linear conductance. 

Fluctuation- dissipation ratio and spin decoherence: A 
number of features underlie the V/T scaling found here 
for the non-linear conductance. First, we observe that 
the kernel in the expression for the current in Eq. ([6]) 
is a function of uj/T and V/T only: — = 

— ^^^rL/iTi ■ This is actually a manifestation of 
the particular form for the fluctuation-dissipation ratio 
of the T-matrix: 



FDR T 



T>+T 



< 



sinh(w/T) 



cosh(w/T) + cosh(V/2T) ' 



(7) 



This result for the fluctuation-dissipation ratio is valid for 
arbitrary bias voltage, regardless of whether the system 
is critical, and for any N. It appears not to have been 
recognized before, but it can be traced back to the steady 
state condition [23]. In light of Eq. (J7J), the V/T scaling 
seen in the I-V characteristics would follow if the spectral 
density of the T-matrix exhibits simultaneous uj/T and 
V/T scaling. The latter is indeed observed. 

To explore the origin of such scaling behavior fur- 
ther, and especially to elucidate the spin decoherence, 
we have studied the local spin susceptibility. We find 
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FIG. 3: Scaling of the local dynamical spin susceptibility. 
Shown are the spectral function (a) and the fluctuation- 



dissipation ratio (b) defined as FDR X 
both taken at uj — V/A. 
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that it too displays a V/T and uj/T scaling. Fig. il- 
lustrates the scaling collapse of the spin spectral density, 
Imx a - It shows the V/T scaling at a particular frequency, 
uj = V/4; independent uj/T and V/T scalings are obeyed 
when the uj and V dependences are separately analyzed 
(not shown). Fig. [SJd shows a similar scaling collapse for 
the fluctuation-dissipation ratio of the spin susceptibility. 
Note that, unlike for T, the steady state condition does 
not fix the fluctuation-dissipation ratio of x- 

Two important points arc in order. First, the vio- 
lation of the fluctuation-dissipation theorem away from 
thermal equilibrium depends on the physical property, as 
had already been noted in other settings (e.g., Rcf. [33j); 
relatedly, the notion of an effective temperature appar- 
ently fails here. Second, the fluctuation-dissipation ratio 
nonetheless still satisfies universal scaling. 

The ui/T and V/T scaling occurs only at the quan- 
tum critical point and in the scaling regime of the 
quantum critical LM phase. As in the equilibrium 
relaxational regime (V = 0, uj -C ksT/h) where 
an ujJT scaling implies a linear- in-T spin relaxation 
rate|2|, in the non-equilibrium relaxational regime here 
(T = 0, uj <C eV/h), the ui/V scaling implies a 
linear-in-U dependence of the dccohcrence rate, Ty = 
H<91nx a (w,T = 0,V)/doj]~l Q = c(e/K)V. It is worth 
stressing that, in contrast to its counterpart in the high- 
voltage (V 3> ksTj^/e) perturbativc regime 

EH EE 13, 

the decohcrence rate here is universal: c is a number 
characterizing the fixed point. 

The Kondo effect involving ferromagnetic leads has al- 
ready been observed and it was argued that imple- 
mentin g th e tuning to the quantum critical regime is 
feasible The scaling of the out-of-equilibrium steady- 
state current can accordingly be studied. Measuring the 
local spin response is more challenging experimentally. 
A promising route could be provided by the single-spin 
ESR technique [35 1 . 

In conclusion, we have identified a system in which 
quantum criticality out of equilibrium can be systemati- 
cally studied. In this model system, we have shown that 
the universal scaling is obeyed by the steady state cur- 



rent, the pertinent spectral densities, and the associated 
fluctuation-dissipation ratios. We have been able to cal- 
culate the entire scaling function for each of these non- 
equilibrium quantities. Our theoretical approach can be 
extended to study the transient behavior and other non- 
equilibrium probes of quantum criticality. 
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Note added- After this paper was posted 
(|arXiv:0805.3717|) . another work (C.-H. Chung et 
al, larXiv:0811.l230| Phys. Rev. Lett. 102, 216803 
(2009)) has appeared on a nonequilibrium dissipativc 
problem with an Ohmic spectrum. Their case has a 
Kosterlitz-Thouless transition, while the sub-Ohmic case 
we study here contains a genuine quantum critical point. 
Their results are complementary to ours. 



[1] Focus issue: Quantum phase transitions, Nat. Phys. 4, 
ppl67-204 (2008). 

[2] S. Sachdev, Quantum Phase Transitions (Cambridge 
Univ. Press, Cambridge, 1999). 

[3] Q. Si et al., Nature 413, 804 (2001); P. Coleman et al., 
J. Phys. Cond. Matt. 13, R723 (2001); H. v. Lohneysen 
et al., Rev. Mod. Phys. 79, 1015 (2007); P. Gegenwart, 
Q. Si, and F. Steglich, Nat. Phys. 4, 186 (2008). 

[4] A. G. Green and S. L. Sondhi, Phys. Rev. Lett. 95, 
267001 (2005). 

[5] D. Dalidovich and P. Phillips, Phys. Rev. Lett. 93, 
027004 (2004). 

[6] K. Damle and S. Sachdev, Phys. Rev. B 56, 8714 (1997). 

[7] A. Mitra et al., Phys. Rev. Lett. 97, 236808 (2006). 

[8] D. E. Feldman, Phys. Rev. Lett. 95, 177201 (2005). 

[9] A. N. Pasupathy et al., Science 306, 86 (2004). 
[10] L. Zhu and Q. Si, Phys. Rev. B 66, 024426 (2002). 
[11] G. Zarand and E. Dernier, Phys.Rev.B 66, 024427 (2002). 
[12] S. Kirchner et al., PNAS 102, 18824 (2005). 
[13] S. Kirchner and Q. Si, Physica B 403, 1189 (2008). 
[14] L. Zhu et al., Phys. Rev. Lett. 93, 267201 (2004). 
[15] N. Wingreen and Y. Meir, Phys. Rev. B 49, 11040 (1994). 
[16] A. Kaminski et al., Phys. Rev. B 62, 8154 (2000). 
[17] P. Coleman et al., Phys. Rev. Lett. 86, 4088 (2001). 
[18] A. Rosch et al, Phys. Rev. Lett. 87, 156802 (2001). 
[19] B. Doyon and N. Andrei, Phys. Rev. B 73, 245326 (2006). 
[20] Y. Utsumi et al., Phys. Rev. B. 71, 245116 (2005). 
[21] J. Martinek et al., Phys. Rev. Lett. 91, 247202 (2003). 
[22] M. Choi et al., Phys. Rev. Lett. 92, 056601 (2004). 
[23] S. Kirchner and Q. Si, in preparation. 
[24] O. Parcollet et al., Phys. Rev. B 58, 3794 (1998). 
[25] D. L. Cox and A. E. Ruckenstein, Phys. Rev. Lett. 71, 
1613 (1993). 

[26] M. H. Hettler et al., Phys. Rev. Lett. 73, 1967 (1994). 
[27] H. Haug and A.-P. Jauho, Quntum Kinetics in Transport 

and Optics of Semiconductors (Springer, 1996). 
[28] A. Kamenev, Lecture notes for the LXXXI Les Houches 



5 



Summer School "Nanoscopic Quantum Transport". 
[29] Y. Meir and N. Wingreen, Phys. Rev. Lett. 68, 2512 
(1992). 

[30] R. M. Potok et al., Nature 446, 167 (2007). 

[31] D. C. Ralph et al., Phys. Rev. Lett. 72, 1064 (1994). 

[32] J. von Delft et al., Ann. Phys. 273, 175 (1999). 



[33] A. Mitra and A. J. Millis, Phys. Rev. B. 72, 121102(R) 
(2005). 

[34] J. Paaske et al., Phys. Rev. B. 70, 155301 (2004). 
[35] M. Kroner et al., Phys. Rev. Lett. 100, 156803 (2008). 



